Spontaneous Symmetry Breaking in Superfluid Helium-4 
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We derive an analytical expression for a critical temperature of spontaneous symmetry breaking in 
a repulsive hard-core interacting Bose system. We show that the critical temperature of spontaneous 
symmetry breaking in a hard-core interacting Bose system is determined by the three physical 
parameters: the density of Bose liquid at absolute zero (po), the mass (m) and the hard sphere 
diameter (a) of a boson. The formula that we have derived is T c = poTvh 2 a/m 2 kB- We report T c 
of liquid helium-4 is 2.194 K, which is significantly close to the A-temperature of 2.1768 K. The 
deviation between the predicted and experimental values of the A-temperature is less than 1%. 

PACS numbers: 67.40.-w 



1. The year 2008 is the 100th anniversary of "Liquid He- 
lium Year" , in which Heike Kamerlingh Onnes produced 
liquid helium. We are in a truly memorial moment in 
the history of low temperature physics. Liquid helium 
has been studied for the past 100 years, and there is 
no question that liquid helium is of central importance 
in this field. (I would like to turn readers' attention to 
Ref. [l[ . It presents the first complete English translation 
of the inaugural speech of Heike Kamerlingh Onnes at the 
University of Leiden in 1882. Although his speech is not 
related to helium, it is quite interesting and gives us a 
lesson. I believe that it attracts the interest of readers.) 

In spite of great efforts, an important problem still re- 
mains to be solved. The problem is as follows. Liquid 
helium-4 undergoes a phase transition, known as the A- 
transition, to a superfluid phase at the A-temperature. 
The experimental value of the A-temperature T\ is ap- 
proximately equal to 2.2 K. The problem is quite simple. 
Why is T\ ~ 2.2 K? What physical parameters determine 
the value of the A-temperature? This is the "Why 2.2 
K?" problem, which is the subject of the present paper. 
It is a long-standing dream in low temperature physics to 
derive a formula for the liquid helium-4 A-temperature. 

In 1938, Fritz London proposed that the A-transition 
of liquid helium-4 probably has to be regarded as the 
Bose-Einstein condensation (BEC) Q. He calculated a 
critical temperature of the BEC in an ideal Bose gas: 



Tbe = 



2nh 2 ( n 
mkB~ U(3/2) 



2/3 



(1) 



where m is a boson mass, n is a number density of Bose 
gas and C(3/2) = 2.612375 • • • . Thus, the Bose-Einstein 
temperature Tbe is determined by the two physical pa- 
rameters, m and n. London reported Tbe of helium is 
3.13 K [3(, which is the same order of magnitude as the 
A-temperature of 2.18 K. His proposal is quite important, 
because it was the first suggestion that the A-transition 
was intrinsically related to the Bose-Einstein condensa- 
tion. Recall that it was not clear in the 1930s whether 
the BEC was an actually physical phenomenon Q . 



Although London's work was revolutionary progress, 
there is a difference of approximately 0.95 K between his 
theoretical prediction and the experimental value. This 
deviation is approximately 44% and it is not negligible. 
It is widely believed that this discrepancy arises from the 
neglect of interactions between helium atoms in his the- 
oretical calculation. This conjecture is quite reasonable 
because liquid helium is not an ideal Bose gas but a re- 
pulsive interacting Bose system. However, in general, it 
is difficult to consider interactions between atoms in an- 
alytical calculations of the BEC. Although this problem 
has been discussed [4|, it remains very important. 

In this paper, we study the "Why 2.2 K?" problem 
via an alternative approach. The key concept in our ap- 
proach is spontaneous symmetry breaking. We derive an 
analytical formula for a critical temperature of sponta- 
neous symmetry breaking in a repulsive hard-core inter- 
acting Bose system. We show that a critical temperature 
is determined by the three physical parameters: the den- 
sity of Bose liquid at absolute zero (po), the mass (m) and 
the hard sphere diameter (er) of a boson. Our formula 
predicts that a critical temperature of liquid helium-4 is 
2.194 K, which is significantly close to the experimental 
value of the A-temperature (T x = 2.1768 K @). 

2. Spontaneous symmetry breaking is a symmetry break- 
ing by the ground state of a system. The symmetry can 
be discrete or continuous, and be local or global. Sponta- 
neous symmetry breaking is a phenomenon that occurs 
in many systems. In condensed matter physics, ferro- 
magnetism is a primary example. For superconductiv- 
ity, the importance of spontaneous symmetry breaking 
was firstly emphasized by Yoichiro Nambu [6J. In parti- 
cle physics, the electroweak symmetry breaking is estab- 
lished in the Standard Model. In cosmology, a number of 
spontaneous symmetry breaking play an important role 
in the history of the early universe. Furthermore, the re- 
cently proposed ghost condensation 0] in the accelerating 
universe can be regarded as a kind of spontaneous sym- 
metry breaking. Thus, spontaneous symmetry breaking 
is a very important concept in modern physics. 
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Spontaneous symmetry breaking is crucial in liquid 
helium-4. In short, the X-transition is the spontaneous 
symmetry breaking. This is the basis of our approach. 
We derive a formula for a critical temperature of spon- 
taneous symmetry breaking in the framework of effective 
field theory described by the order parameters. In the 
case of liquid helium-4, the order parameter is a macro- 
scopic wavefunction ip, which is a one-component com- 
plex scalar field with both amplitude and phase. We 
begin with the theory described by this scalar field <p. 
3. The Lagrangian density £(= K — V) of the system is 
given by the nonrelativistic Goldstone model of the form 



V 



A 



-fnp ip 



A 



2nh 2 a 



m 



(2) 



The kinetic part JC = ihp*dtp + h 2 p*V 2 p/2m is not 
important here, and therefore we concentrate on the po- 
tential V. In Eq. p is the chemical potential, A is the 
coupling constant for helium interactions and ip* is the 
complex conjugate of ip. The interaction between helium 
atoms is repulsive at short distances (hard-core interac- 
tions). Therefore, the coupling constant A = 2irH 2 o-/m 
is positive, where to and a are the mass and the hard 
sphere diameter of a helium atom, respectively. 

We first consider the physical dimensions of parame- 
ters of the potential In units of ft = 1 and fcs = 1 
(ft is the reduced Planck's constant and fc# is the Boltz- 
mann constant), for any physical quantity Q, its physical 
dimension [Q] is written as the product of the dimensions 
of length and temperature: 



[Q] = [Length]" [Temperature] 



(3) 



where a and (3 are numerical constants. The physical di- 
mension of the Lagrangian density is [£] = [L] ~ 3 [T] , and 
that of the scalar field <p is [up] = [L]~ 3 / 2 . (Here [L] and 
[T] represent [Length] and [Temperature], respectively.) 
Hence it is straightforward to derive the physical dimen- 
sions of parameters in Eq. 

M = [T], [A] = [L] 3 [T], [to] = [L]- 2 [T]-\ [a] = [L]. (4) 

Note here that in the present case the coupling constant 
A is a dimensionful parameter, in contrast to a relativistic 
case in which A is dimensionless. These physical dimen- 
sions help us to understand the physics of liquid helium-4. 

We next consider the symmetry of Lagrangian. The 
potential V of Eq. @ is invariant under the U(l) phase 
transformation of the field p: p — > e l6 (p, p* — ► e~ %e p* . 
Furthermore, if this is a global transformation in which 
the transformation parameter 9 does not depend on the 
space and time coordinates, the kinetic part K is in- 
variant under the U(l) phase transformation. Hence, 
the system described by the Lagrangian density C has 
a global U(l)~0(2) symmetry. This global U(l) symme- 
try is spontaneously broken at low temperature as you 
see below. 
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FIG. 1: The potential V in the case that fj, > 0. 

4. The ground state of a system is given by solving the 
condition dV /dp — 0. First, we consider the case that 
a chemical potential p is negative (jj, < 0). In this case, 
there is only a trivial solution p = 0. This solution is 
the stable ground state, because d 2 V /dp*dp\ ip= o > 0. 
Hence the global U(l) symmetry is unbroken in this case. 
Second, we consider the case that p is positive (p > 0). 
In this case, in contrast to the case that p < 0, there is a 
nontrivial solution in addition to a trivial solution ip = 0: 



(5) 



m = \ t 



Fig-fflshows the potential V in the case that p > 0. In the 
present case, the trivial solution p> = is unstable and 
it is not the ground state of the system. The nontrivial 
solution \p\ = \f n/\ is the stable ground state, because 
the curvature of the potential at two solutions satisfies 



d 2 v 



dp*dp 



tp=0 



<0, 



d 2 V 



dp*dp 



> 0. 



(6) 



Therefore, if a chemical potential p is positive {p > 0), 
the global U(l) symmetry is spontaneously broken. From 
above arguments, we reach a following view. A physical 
system undergoes a phase transition at a critical temper- 
ature T c , if a chemical potential satisfies the conditions: 
p < (T > T c ), p > (T < T c ). Spontaneous symme- 
try breaking occurs only if a chemical potential becomes 
positive. Therefore, in contrast to the Bose-Einstein con- 
densation, spontaneous symmetry breaking never occurs 
in an ideal Bose gas, in which a chemical potential is 
always negative (p < 0). This is a crucial difference be- 
tween spontaneous symmetry breaking and the standard 
Bose-Einstein condensation. 

5. Before deriving a formula for a critical temperature, 
we consider the temperature dependence of a chemical 
potential p. Although the temperature dependence of p 
is not necessary for deriving a formula, it is useful to un- 
derstand a phase transition. For this reason, we consider 
the temperature dependence of p in the following. The 
quantity \p\ 2 is equivalent to the number density n of 
superfluid. Therefore, the following relation is satisfied 
below the A-temperature 



m 




(7) 



where p s is the superfhrid density and mis a boson mass. 
From this relation, the superfhrid density is given by 



Ps 



mp 



(8) 



This expression indicates that the temperature depen- 
dence of p is obtained from that of p s . 
Here we take a minimal model: 



= 1 - 



(9) 



where p{— p n + p s ) is the total density of liquid helium, 
and p n , p s are the normal fluid and superfhrid densities, 
respectively. It should be emphasized that this model 
does not affect a formula for a critical temperature. Fig. [2] 
shows the superfiuid density ratio as a function of tem- 
perature. Two functions 1 - (T/T\) 5 and 1 - (T/T\) 7 
are plotted, for comparison. Fig. [2] shows that Eq. (J9j) 
and experimental values are consistent except near the 
A-transition, at which the superfiuid density ratio is of 
the form (1 - T/T A ) 2 / 3 . The form of Eq. (9) is the same 
as that of the Bose-Einstein condensation in an ideal Bose 
gas, 1 - (T/T c ) 3 / 2 . Although the exponent "6" has not 
yet been derived, we do not consider this problem here. 
From Eqs. ((8]) and (|9|), we obtain the expression 



M = Mo 1 - 




Pa = 



(10) 



where po is a chemical potential at absolute zero and po is 
the total density at absolute zero. We have used the ap- 
proximation p ~ pq in Eq. (|10[) . because the total density 
p is approximately a constant below the A-temperature. 
Therefore, Eq. (fTTj|) is a good approximation at low tem- 
perature (T < 1.7 K). Fig. [3] shows the schematic tem- 
perature dependence of the potential V. For T > T\, the 
potential has only one minimum at ip — and the curva- 
ture at ip = is positive. When T — T\, the curvature at 
the minimum is zero, d 2 V /dp*dip\ v= o = 0. For T < T\, 
the curvature at ip — is negative and the scalar field 
rolls down to the \(p\ ^ minimum. Thus, there is no 
barrier of potential between the minimum at ip = and 
l^l =/= 0, indicating a second-order phase transition. 
6. We now derive a formula for a critical temperature of 
spontaneous symmetry breaking. When T — 0, the field 
ip is located at the nonzero minimum |iys| 0. Hence the 
global U(l) symmetry is broken in this phase. As the 
temperature increases above a critical temperature, the 
potential has only one minimum at ip — and the global 
U(l) symmetry is recovered. From these observations, we 
reach a following conclusion; the depth of the potential 
well at absolute zero determines a critical temperature. 
When T = 0, the depth of the potential well is given by 



v(M = vWA) 



p n 



(11) 
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FIG. 2: Superfiuid density ratio for liquid helium-4 as a func- 
tion of temperature. The inset shows the temperature region 
1.0 K < T < 1.7 K. Key: (dotted line) 1 - (T/T A ) 5 ; (solid 
line) 1 - (T/Txf; (dashed line) 1 - (T/T A ) 7 ; (+) the data 
from Ref. d. 




FIG. 3: The schematic temperature dependence of the poten- 
tial V for a second-order phase transition. 



where n$ — po/m is the number density of Bose liquid 
at absolute zero. This quantity pquq/2 represents the re- 
quired energy density to recover the global U(l) symme- 
try of the system. Therefore, the corresponding thermal 
energy density fcsT c n,o is equivalent to pqUq/2. Conse- 
quently, we obtain the following relation 



knTc 



Mo 
2 ' 



(12) 



Fig. 0] shows the present situation. From above argu- 
ments, we have the following relations: 



, rr, fJ-o Po Po < 
2mA 



2irh 2 a 



(13) 



Therefore, a critical temperature of spontaneous symme- 
try breaking in a repulsive interacting Bose system is: 



poirh 2 o- 
m 2 kB 



(14) 



where po is the density of Bose liquid at absolute zero, 
m and a are the mass and the hard sphere diame- 
ter of a boson, respectively. Substituting the values of 
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FIG. 4: The absolute zero potential per number density. 

liquid helium-4 into Eq. [fl4")). p = 0.1451 g/cm 3 @, 
m = 6.6465 x 10~ 24 g, and a = 2.639 A [1,0, we obtain 
the value T c = 2.194 K. This prediction is significantly 
close to the experimental value of the A-temperature 
2.1768 K [f|. The deviation A between the predicted and 
experimental values of the A-temperature is less than 1%: 



A = 100 x 



T, 



(theory) 



r. 



(exp) 



T- 



(exp) 



< i%. 



(15) 



Thus, our prediction is much closer to the experimen- 
tal value than London's prediction. It is because in the 
present approach the repulsive interactions between he- 
lium atoms are included in the derivation of a formula. 
While the Bose-Einstein temperature Tbe in an ideal 
Bose gas is determined by the two physical parameters, 
our formula for a critical temperature, Eq. (|14p . is deter- 
mined by the three physical parameters: the density of 
superfluid (po), the mass (m) and the hard sphere diam- 
eter (<r) of a boson. The additional physical parameter a 
includes information about interactions, and it has signif- 
icantly improved the theoretical prediction for a critical 
temperature. The results are summarized in Table HI 

Finally, we comment on the hard sphere diameter a. 
Although po and m have been determined precisely, the 
uncertainty of a is relatively large. The hard sphere di- 
ameter a represents the point at which the interatomic 
potential V(r) is zero (V(a) = 0), and the distance r m 
represents the point at which the potential V(r) is min- 
imum (dV (r) I dr\ r=rm = 0). Although the value of r m 
has been reported by many groups [l(| , that of a has been 
less reported. Furthermore, the relation r m = 2 1 / 6 ct of 
the (12,6) Lennard-Jones potential model is not appli- 



TABLE I: A comparison between BEC and SSB. Right row 
is given in this work. The deviation A is defined by Eq. (|15J| . 




1.58 2.6 2.62 2.64 2.66 

the hard sphere diameter of a helium atom [A] 

FIG. 5: The A-temperature as a function of the hard sphere 
diameter a of a helium atom. 



cable to the present study, because it is not sufficiently 
a good approximation. For these reasons, we plot the 
A-temperature as a function of a. Fig. [5] shows that the 
deviation between the predicted and experimental values 
is less than ±1% for the region 2.592 A < a < 2.644 A. 
7. In conclusion, we have derived an analytical for- 
mula for a critical temperature of spontaneous symmetry 
breaking in a repulsive interacting Bose system. The for- 
mula that we have derived, Eq. (|14[) . is a simple analytical 
expression to predict T\ ~ 2.2 K. We hope that this work 
contributes to progress of low temperature physics. 



System 


ideal Bose gas 


interacting Bose liquid 


[9] 


Key concept 


BEC 


SSB 




T c 


27r/j 2 (n/C(3/2)) 2/3 /fnfcs 


poTrh 2 a/m 2 kB 


[10] 


parameter^ 


2 


3 


T c ( 4 He) 


3.13 K 


2.194 K 




A 


44 % 


0.8 % 
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